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Abstract 

An exact analytic form for the second virial coefficient, valid for the entire range of temperature, is 
presented for the Lennard-Jones fluid in this paper. It is derived by making variable transformation 
that gives rise to the Hamiltonian mimicking a harmonic oscillator-like dynamics. It is given in 
terms of parabolic cylinder functions or confluent hypergeometric functions. Exact limiting laws 
for the second virial coefficient in the limits of T — > and T — > are also deduced for the 
Lennard-Jones fluid. They have the forms: — 16\/27rt)oe^^ as T — > and 4-v/2r (|) vq 

as T — > oo, where e is the well depth and vq = ira^ /Q with a denoting the size parameter of the 
potential, and j3 = l/ksT. 
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I. INTRODUCTION 



The statistical mechanical formula for the second virial coefficient can be easily com- 
puted by using a numerical method or, if temperature is sufficiently high, by a series ex- 
pansion method^. As a matter of fact, a numerical table is available^ for the Lennard- Jones 
fluid. Therefore it does not pose a practical problem, although a low temperature expan- 
sion method is not available. Nevertheless, it would be interesting from the theoretical and 
pedagogical as well as aesthetic standpoints and also for the practical utility, if there were 
available exact analytic results for realistic interaction potential models, which are valid for 
the entire range of temperature. In the Iiterature2'i2., an exact analytic form for the second 
virial coefficient for the Lennard- Jones (LJ) fluid was obtained by using coordinate trans- 
formation in the cluster integral and reading off the integral tabled to obtain such a form. 
However, such transformed integrals have a deeper underlying dynamics that mimics a dy- 
namical system obeying a harmonic oscillator potential, but reading off an integral table 
does not reveal the underlying dynamical structure therein. Therefore, if such a feature is 
made evident, one can gain a considerable insight into the dynamics of the LJ fluici^. The 
method employed also provides a valuable lesson on how to handle such integrals that might 
appear in the study of statistical mechanics of simple liquids. This method is not available 
elsewhere in the literature as far as this author is aware of. 

In this paper, we present an exact analytic result for the second virial coefficient of the 
Lennard- Jones (LJ) fluid, which is obtained without using an expansion method and valid 
for the entire range of temperature. Exact limiting forms are also deduced therefrom as 
T — > and T —>■ oo. The second virial coefficient obtained is given in terms of parabolic 
cylinder functions or confluent hypergeometric functions, which are convergent and well 
deflned for all values of temperature. The form presented for the second virial coefficient 
therefore is valid for all temperatures. 

II. ANALYTIC SOLUTIONS 

The second virial coefficient^ of the LJ fluid may be written in the reduced form 
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where vo = na^/G, the volume of the contact sphere of diameter a {a = size parameter of 
the LJ potential), e is the well depth, and (3 = l/ksT, inverse temperature, being the 
Boltzmann constant. The object of interest is the integral / in Eq. ([T]). With transformation 
of variables 

a = ^/ep, (2) 

y = — ' (3) 



the integral / can be put into the form 

na) = ^Jia), (4) 
where J (a) is defined by the integral 

POO 

J{a)= / dyy-'/^ (e-^^e2"^/^ - l) . (5) 
Jo 

This integral is usually evaluated either by a series expansion method or by a numerical 
method. Luckily, this integral is in a form readily available from the integral table^ in a 
known functional form. However, it can be analytically evaluated without using an integral 
table or expansion method, as will be shown below. 

To achieve this aim, perform integration by parts once to obtain J (a) in the form 

J (a) = 4 [aBy^ (a) - B^/^ (a)] , (6) 

where i?i/4 and i?3/4 are defined by the integrals 

POO 

i?i/4 (a) = / dyy-^/^e-ye^'^^, (7) 
Jo 

POO 

53/4 («) = / dyy-'/^e-ye^''^. (8) 







These integrals are functions of parameter a. Differentiating these integrals with a, we 
obtain a pair of first-order differential equations 
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da 2 



^3/4 («) + 2a5i/4 («) , (9) 



dB 



- 2Si/4 («) . (10) 



da 

This pair of differential equations can be combined to a single homogeneous second-order 
differential equation: 

d'^ih dib , . 



with the simphfied notation 



With the transformations 



and 



(a) = 5i/4 (a) . 
z = \/2a 



ei^'^d (z) 



the differential equation (fTTl) can be transformed into a standard form 
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(z) = 0. 



(12) 
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This is akin to the Schrodinger equation for a harmonic oscillator — of a negative eigenvalue 
in the present case. Therefore, it represents a particle of negative energy subjected to a 
parabolic potential. It, in fact, is a differential equation for parabolic cylinder functions^i^. 
Its two independent solutions, one even and the other odd function of z, are: 
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01 (z)=e-t^^M(-,-,-z2 
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Here M(a, 6, t) is a confluent hypergeometric function of Kummer- 



M{a,b,t) = Y, 



n=Q 



(6)„ n\ 



(18) 



where 



(a)o = 1' 

(a)„ = a (a + 1) (a + 2) ■ ■ ■ (a + n - 1) (n > 1) 



(19) 



It is convergent for all values of t. Its asymptotic form will be interest to us later: for positive 
real t it is given by 

m—l / 



M(a,6,t) = ^e*t"-^ 



n=0 



(20) 



This formula may be used to compute the solutions for a large value of t at fixed values 
of a and h. The solutions 0i {z) and 02 (-z), in fact, are parabolic cylinder functions, which 
are generic solutions for the Schrodinger equations for quadratic potentials. This implies 
that the dynamics of the LJ potential fluid closely resembles that of a harmonic (quadratic) 
potential. 

Therefore the general solution for ip{z) is: 

1 1 2\ 3 1 

, -, -z^ + C2zM -, -, 
^4' 2' 2 ) ^ V4 2' 2 

where Ci and C2 are constants, which may be determined by considering the boundary 
conditions. 
Noting that 



Bi/i{z) = ij{z) = ciM 



(21) 



-f-M(a, b, t) = yM{a + 1, 6 + 1, t), 
at b 



(22) 



we find 



J{a) = 4ci 

+ AV2c2a 
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The coefficients ci and C2 can be determined by examining the limiting form of J {a) as 
a (a boundary condition). From Eq. fl23|) 



J(a) 



-ci 



+ y2c2a + O (a^) 



(24) 



whereas direct evaluation of J{a) by series expansion of the factor exp {2a^^ in Eq. ([5]) 
yields 

J(«) = -4r (^^) + 2r (i) a + O [a') , (25) 

where F (1) and F (|) are gamma functions: F (1) = 3.62560 ■ ■ ■ and F (|) = 1.22541 • • • . 
Comparing Eqs. fl2^ and fl25l) . we find 

''^ (26) 

(27) 



ci = r(- 



C2^^Ffi,. 



2v/2 V4 



Thus J (a) is now fully determined. 
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Putting together the results produced up to this point, we finally obtain the second virial 
coefficient in the form 
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This is the result we have set out to show for the LJ fiuid. One may try to put this result 
into a simpler form by using the recurrence relations of Kummer's functions, but the present 
form appears to be an optimum form. Rigorous limiting laws can be deduced for B2 from 
Eq. m- 

The limiting form of -B2 as T — > cxo or — is easily deduced to be 

'3^ 



Bo 



4v^r - Uo {epf [l + 0{eP)]. 



(29) 



Thus B2 +0 as T ^ 00. This means that there is a high temperature regime where B2 
is positive, and as T ^ 00, it vanishes on the positive side according to the law indicated. 

On the other hand, the limiting form of B2 as T ^ or 6j3 —>■ 00 is deduced from the 
asymptotic forms of the confiuent hypergeometric functions given in Eq. (!20l) . We find 



B2 (T) = -16V27rwoe"^ {ep)^ 
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(30) 
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This limiting law shows that B2 (T) tends to negative infinity according to the formula 
indicated and is negative below a certain point in T. 

These limiting laws for B2 are not easily deducible from Eq. ([l]) or Eq. ([5]) or the 
series expansion form^ thereof, but they are simple to deduce if the exact analytic solution 
presented is made use of. 

From the limiting behaviors fl29|) and fl30|) we can conclude there must exist a point in 
T at which i?2 (T) crosses the T axis (i.e., becomes zero), that is, the Boyle temperature 
Tb = PB^/ks is defined, as usual, by 

B2 (Tb) = 0. (31) 

According to the analytic result obtained, the Boyle point is determined from a real root of 
the equation 
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Its numerical solution yields 

T* = (epB)'^ = 3.41793, (33) 

which should be compared with the literature value^ = 3.42. This value is practically 
attained with the truncation of M (a, b, e^b) at n = 3. 

In conclusion, we have presented an exact analytic second virial coefficient of the LJ fluids, 
which are valid for the entire temperature range, and its asymptotic behaviors (limiting laws) 
as T — > or T oo. In view of the agreement of the Boyle temperature with the literature 
value deduced from the table for the second virial coefficient^ it seems to be unnecessary to 
tabulate the numerical values of the second virial coefficients; it is rather trivial to do so. 
The utility of the result obtained is self-evident for some deductions one can make about 
thermodynamic properties of the LJ fluid. Compared to the method that simply reads off 
the integral table upon variable transformation in the integral for B2, the present method 
provides considerable insights into the dynamics of the LJ liquid. 

Finally, it is useful to note that the present exact analytic result for the second virial co- 
efficient owes its existence to the mathematically favorable combination of the exponents 12 
and 6 of the potential that produces the closed form for the differential equation for (a), 
Eq. (ITT!) . For other potential models consisting of repulsive and attractive branches with 
different exponents we do not obtain a closed differential equation, but a open hierarchy of 
ffist-order differential equations for integrals making up J (a). The case of exponents (9, 6), 
namely, the L J (9, 6) potential, produces a closed inhomogeneous second order differential 
equation, but its solutions do not seem to be simple and clean. 
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